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Abstract. I review the general aspects of cosmological parameter estimation from 
observations of the cosmic microwave background (CMB) temperature anisotropics in 
the framework of inflationary adiabatic models. The most recent CMB datasets are 
starting to give good constraints on the relevant parameters of inflationary adiabatic 
models. They point toward a model consistent with the basic predictions of inflation: 
a nearly flat universe, with a nearly scale invariant spectrum of primordial fluctuation. 



INTRODUCTION 

The study of the cosmic microwave background (CMB) temperature anisotropics 
has long been recognized as one of the most powerful tools to answer the basic ques- 
tions about the nature of the universe: what is its geometry, its matter and energy 
content, what are the initial conditions which seeded the formation of structure, 
etc. It is a firm theoretical conclusion that the angular distribution of the CMB 
anisotropies must encode a vast amount of information on the cosmological param- 
eters. The majority of this information is thought to be concentrated at angular 
scales smaller than about 1 degree on the sky, corresponding to regions of the uni- 
verse that were in casual contact when the background photons decoupled from the 
matter (at redshifts of about 1000). On this scales, physical processes in the early 
universe were able to leave an imprint on the CMB. 

The pioneering observations by the COBE satellite in the early 90's, which led 
to the first unambiguous detection of the CMB anisotropies [1], were followed by a 
large number of ground based and balloon-borne observations which attempted to 
collect information on the fine-structure pattern of the anisotropy, where most of the 
dependence on cosmological parameters is encoded. Recently, the BOOMERanG 
[2] and MAXIMA [3] balloon-borne experiments produced the first high-resolution 
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maps of the CMB anisotropy pattern, and moved us closer to the goal of a long- 
awaited high precision measurement of the CMB angular power spectrum. These 
results have been used by several authors [4] [5] [6] [7] to obtain high-precision 
constraints on the set of cosmological parameters which defines the inflationary 
adiabatic class of models. 

In this review I will attempt to give an idea of how the cosmological parameters 
affect the CMB observables, and of the process which leads from the observation of 
the CMB anisotropy to the extraction of such parameters. The current constraints 
on inflationary models from the CMB will also be discussed. 

THE CMB ANISOTROPY DEPENDENCE ON 
COSMOLOGICAL PARAMETERS 

It is generally believed that the observed large scale structure of the Universe 
formed by gravitational amplification of small density perturbations generated 
in the early universe. In such gravitational instability scenarios, the presence of 
anisotropies in the temperature distribution of the CMB is unavoidable: density 
fluctuations must leave an imprint in the CMB at the time of photon-matter de- 
coupling, at redshifts of about 1000. 

The anisotropy as a function of the direction of observations can be expanded in 
spherical harmonics: 

^(7) = E°^WM). (i) 

The coefficients Ce = (\ae m \ 2 ) define the angular power spectrum of the CMB 
anisotropy. The C/s do not depend on the azimuthal index m as a consequence of 
the isotropy of space. For Gaussian initial conditions, the angular power spectrum 
Cg carries all the information on the angular temperature anisotropy of the CMB. 
Each £ probes an angular scale 9 on the sky given approximately by I ~ tt/0. Only 
ts corresponding to angular scales which were in casual contact at decoupling 
may have been affected by physical processes prior to decoupling. For this reason 
the dependence on physical parameters is mostly found at high ts (small angular 
scales), while low ts probe the primordial shape of the power spectrum 2 . 

Within the inflationary adiabatic family, a given cosmological model is specified 
by the value of a number of parameters. These include the fractional density of 
matter in the universe, which is the sum of contributions from baryons and cold 
dark matter, Q m = Qb+Q c d m ; the fractional density of vacuum energy, the total 
energy density, Q = Q m + Q\, which defines the curvature of the universe through 
Qk = 1 — Q; the Hubble constant, parameterized by its value h in units of 100 
km s _1 Mpc -1 ; and the amplitude A and spectral index n of the primordial power 
spectrum of density fluctuations, modeled as Ak n (but see [8] for an alternative 



2 ) Neglecting secondary processes which may alter the spectrum after decoupling. 



view). The C/s corresponding to a set of parameters can be computed exactly 
using high-accuracy numerical codes [9]. A simplified analytical treatment was 
used by some authors (see e.g. [10] for an excellent review on the subject) in order 
to give a better intuition of how different physical processes leave an imprint on 
the CMB angular power spectrum. 

Before decoupling, the photons and baryons are tightly coupled by different scat- 
tering processes. The cold dark matter is non-interacting and contributes only to 
the gravitational potential. The dynamics of the photon-baryon fluid is described 
by an equation reminiscent of the classical Jeans equation governing the perturba- 
tion in a self-gravitating gas, which for a given wave-number k is: 

where a is the scale factor describing the expansion of the universe, R = 3p#/4p 7 
is the baryons to photons density ratio and c s = c/[3(l + R)] 1 ^ 2 is the sound ve- 
locity. For adiabatic initial conditions (i.e. conserving the entropy of the radiation 
per baryon) the intrinsic temperature fluctuation of the photons, ST/T is related 
to the matter density perturbation simply by 5T/T = 5/3. The dots represent 
derivatives with respect to the conformal time r] = / dt/a. The term F describes 
the gravitational effects and can be held approximately constant near decoupling. 
So, the evolution of the perturbations prior to decoupling is essentially governed 
by a forced harmonic oscillator equation. The expansion of the universe introduces 
a viscosity term through a/a that may be neglected for the sake of simplicity. 

The physical interpretation of this equation is very simple. The baryons tend to 
collapse due to self gravitation. The restoring force is provided by the radiation 
pressure k 2 c 2 /3. This sets up acoustic oscillations (the sound velocity c s quantifies 
the resistance of the fluid to compression). The higher R, the larger the amplitude 
of the oscillations. The driving force term due to gravitation, constant in our 
approximation, simply displaces the zero point of the oscillations. Increasing R 
(i.e. the baryon content of the universe) enhances this displacement, and gives more 
amplitude to compressions over rarefactions, because of the increased inertia of the 
fluid. If we freeze the oscillations at the time of decoupling r] dec , each mode will be 
in a different stage of oscillation. The total power will have the largest contributions 
from modes having kc s f]d ec — kr^ ec = mir, where r^ ec = c s i]d ec defines the physical 
scale of the sound horizon at decoupling. This results in a harmonic series of peaks 
in the angular power spectrum, whose position is related to 9 S , the characteristic 
angular scale subtended by the sound horizon at decoupling, by I ~ ir/0 s . Odd 
peaks are due to compression of the fluid, even peaks to rarefaction: so the odd 
peaks will be generally higher than the even peaks because of R. Increasing the 
baryon content will enhance this effect. 

To calculate the dependence of the peaks position on cosmological parameters we 
have to specify the angular diameter distance relation which maps a given physical 
scale at decoupling into an angle 9. As we saw, the relevant physical scale is the 



sound horizon at decoupling rf ec = c s r] dec . The angular scale of the sound horizon 
at decoupling, 9 S , is approximately given by: 

sinhJ fn m \ 1/2 f 1 c nfda 

S ~ 2c s a l H U J ' Ja dec [a m a + n A a* + Q fc a 2 ] V* ' { ) 

for an open universe, flk > 0; 

x sin/ /^y/ 2 yi c (-^Q 1 / 2 ^ 

S ~ 2c s alH ' [fimO + Q A a 4 + 1W] 1 / 2 ' ^ > 

for a closed universe, flk < 0, and: 

- ^ ( ft »,)" 2 ; / = /' r iL #nTT7v (5) 



for a flat universe, = 0. 3 

This very simplified discussion gives an intuitive idea of how the main CMB 
observables, namely the position and height of peaks in the Cgs, are affected by 
cosmological parameters in inflationary adiabatic models. For a given primordial 
power spectrum, the anisotropy pattern is defined by the baryon content and dark 
matter content, Vt^h 2 and Q c d m h 2 , which basically affect the amplitude of fluctu- 
ations at different physical scales and fix the relative heights of the peaks. This 
structure is mapped into different angular scales depending on a combination of 
Q m , £l\ and f2 fc , which affect the position of the peaks. These effects are illustrated 
in Figure 1. This also shows that the parameters actually enter in the power spec- 
trum in a combined way, leading in some cases to almost exact degeneracies: a 
dramatic example of this behavior is shown in Figure 2. The problem of parameter 
degeneracies has been explored in great detail in [11]. 



COSMOLOGICAL PARAMETER ESTIMATION FROM 

CMB MEASUREMENTS 

As we saw in the previous Section, a given cosmological model is specified by a 
set of cosmological parameters: 

P = {PuP2,---,Pn} (6) 

entering in the calculation of the theoretical power spectrum, Cg. A CMB experi- 
ment measures the temperature fluctuation of the CMB in different directions on 
the sky. These data are used to build a minimum variance map of the CMB. From 
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FIGURE 1. The angular power spectrum of the CMB for inflationary adiabatic models with 
scale invariant primordial fluctuations (n=l). The three curves show the effect of varying the 
value of the parameters p= (£1 = f2 m + Ct\, flbh 2 , fl c d m h 2 , Oa) is shown in the three curves. The 
solid curve has p=(l, 0.03, 0.2 0.54). If fif,/i 2 and Q c d m h 2 are kept fixed while varying Q m and 
Oa, the peaks are simply shifted (dotted curve, with p=(0.46, 0.03, 0.2, 0)). When we fix Cl m and 
Q\ and we vary the baryon and cold dark matter content, the positions are basically unchanged, 
and the ratio of peaks height changes (dashed curve, with p=(l, 0.025, 0.17, 0.54)). The points 
are the bandpower measurements from BOOMERanG and MAXIMA, re-scaled within their la- 
calibration uncertainty. 

the map, maximum likelihood estimates of the power spectrum are extracted as a 
set of average bandpower measurements over an interval b in £: 



where the shape function Sg is usually assumed to be Se oc £(£ +1). 

The best estimate of the cosmological parameters from a set of measured band- 
powers can be derived in a Bayesian sense by maximizing the likelihood function: 



We should note that the likelihood is not a Gaussian function of the C b 's. One way 
to estimate how the likelihood depends on the C^'s is to use the ansatz described 
in [12] where an additional quantity x b related to the noise of the experiment is 
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FIGURE 2. Degeneracy lines in the Q m -£lA plane, assuming h = 0.7 and fife/i 2 = 0.03. The 
lines correspond to models having approximately the same angular scale of sound horizon at 
decoupling, and are labeled by the approximate position of the first peak in the CMB angular 
power spectrum. The ticker lines give a reference interval corresponding to 180 < l pea k < 220. 
Flat models lie on the diagonal dotted line (fi = 1). The shaded curve represents the 95 % 
confidence levels from high redshift type la supernovae. 



needed in order to fully characterize the likelihood. This has been shown to work 
well when compared to a brute force exact calculation of the likelihood. 

In a Bayesian framework, getting the a posteriori likelihood of one (or more) 
parameters out of the full likelihood, involves an integration (marginalization) of 
the likelihood over the unwanted (nuisance) parameters. The full likelihood has to 
be weighted by the a priori likelihoods of the nuisance parameters: 

£(P*» • • • » Pi) = f dpi--- dp m C(pi, •••,Pj\pi,---, Pm)C(pi) ■ ■ ■ £(p m ) (9) 

This complicates the analysis considerably, making the problem strongly non-local. 
The overall shape of the likelihood influences the result on subsets of parameters, 
and we then need to characterize the likelihood over the entire parameter space 



TABLE 1. Current constraints (68% c.l.) from CMB data only. 



Dataset ft fl^h 2 n c< i m h 2 n 

BOOMERanG + COBE 1.15g;Jg 0.036g;j$f 0.24g;g| 1.04g; 

MAXIMA + COBE 1.00g;?| 0.030g;8g| 0.20g;^ 1.08g : : 

MAXIMA + BOOMERanG + COBE l.ll&g? 0.032°;°°;! 0.14g;gg l.Olg; 



even if we are just interested in a subset of parameters. The problem is exacer- 
bated by the presence of correlations in the parameter space. In addition, we need 
a knowledge (or a reasonable guess) of the a priori likelihood of the nuisance pa- 
rameters. This can be chosen to be uniform (in the absence of previous knowledge) 
or can be derived from other observations. We also stress the fact that the extent 
of the parameter space effectively acts as a top-hat prior, cutting out a priori some 
values of the parameters. In general one should carefully explore how changing the 
priors on some parameters affects the results for the others. This is particularly 
important when the results for some parameters conflict with those coming from 
other observations. The issues related to the use of different priors have been thor- 
oughly explored in the analysis of the latest CMB datasets (see, e.g. [4] [5] and 
[7])- 



CURRENT CONSTRAINTS 

Over the past decade, many experiments have gathered data on CMB anisotropy, 
which, taken collectively, gave an indication of the shape of the angular power spec- 
trum. Such data were used by a number of authors to set constraints on the param- 
eters of the inflationary adiabatic model (see, e.g. [13], and references therein). The 
quality of CMB data has considerably improved after the recent BOOMERanG 4 
[2] and MAXIMA 5 [3] balloon-borne missions. BOOMERanG estimated the power 
spectrum in 12 bins, over the range 26 < £ < 625, from a map of a 1800 square de- 
grees patch of the southern sky; MAXIMA estimated the spectrum in 10 bins, over 
the range 36 < £ < 785, from a map of a 124 square degrees patch of the northern 
sky. The two measurements are in remarkable agreement and show unambiguously 
the presence of a sharp peak in the region 180 < £ < 220 (see Figure 1). This 
is, in itself, a strong evidence in favor of inflationary adiabatic models; alternative 
theories either predict a broader peak at higher £ or a broad shelf at £ < 200 (see 
e.g. [14]). 

The BOOMERanG and MAXIMA datasets (in combination with COBE) were 
used independently to set constraints on a seven-dimensional space of cosmological 
parameters [4] [5]; a joint analysis was also carried on [7]. Using no prior information 
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TABLE 2. Current constraints (68% c.l.) from CMB data plus prior 
information from high redshift type la supernovae 



Dataset 


n A 


£7 m 


BOOMERanG + COBE + SNe la 

MAXIMA + COBE + SNe la 

MAXIMA + BOOMERanG + COBE + SNe la 


o.728;8I 
0.608:8? 
0.758:8? 


u -°'0.07 
u -°'0.07 
u - oo 0.07 



(except the constraint that the universe is older than 10 Gyr and that the Hubble 
constant is 0.45 < h < 0.90) the analyses performed in [4] [5] and [7] found the 
results reported in Table 1. These CMB data alone are already giving a consistent 
picture, which is in agreement with the basic predictions of inflation: the universe 
is nearly flat, and the primordial fluctuations have a nearly scale-invariant power 
spectrum. Moreover, there is indication of a substantial contribution from non- 
baryonic matter. It may also be observed that the best estimate of the baryon 
density from CMB shows some tension with the value from big bang nucleosynthesis 
[15] Vt b h 2 = 0.019±0.002. Some implications of this discrepancy have been explored 
in [16]. 

As we saw in the first Section, models with different values of fl m and fl\ may 
result in the same angular scale of sound horizon at decoupling. As a consequence, 
it is hard for CMB data alone to separately determine fl m and fl\. This is clear 
from Figure 2 where we plot lines in the Q m — Q\ plane corresponding to the same 
angular scale of sound horizon at decoupling. To break this degeneracy, prior 
information from different datasets can be used. The results obtained by [4] [5] and 
[7] when combining the CMB data with constraints in the Q m — Q\ plane coming 
from observations of high- redshift type la supernovae [17] [18] are shown in Table 2. 
These results make a very strong case for the existence of a substantial contribution 
from some form of unknown negative-pressure component (named "dark energy" in 
the recent literature). Attempts to investigate the nature of this component using 
the CMB, while assuming strong priors on the other parameters, have been made 
in [19]. 

Finally, it is remarkable that, as shown in [7], when completely independent priors 
from large scale structure observations are used instead of those from supernovae, 
totally consistent results are obtained. 

The CMB is already proving very powerful in improving our knowledge of cos- 
mological parameters. Increasingly accurate measurements of the power spectrum 
will come in the next few years by satellite missions such as MAP 6 and Planck 7 , 
which will further strengthen our understanding of the nature of the universe. 
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